home *** CD-ROM | disk | FTP | other *** search
/ The Complete Utilities To…ka 501 Killer Utilities! / 501 Killer Utilities! (Macworld July 1995).cdr / Programming / RLaB / examples / sfft.r < prev    next >
Encoding:
Text File  |  1994-09-20  |  355 b   |  26 lines  |  [TEXT/RLAB]
  1. //
  2. // Calculate Fourier transform of a finite signal (SLOWLY)
  3. //
  4.  
  5. sfft = function ( x )
  6. {
  7.   global (pi)
  8.  
  9.   if(class (x) != "num") 
  10.   {
  11.     error ("input to sfft must be a matrix");
  12.   }
  13.  
  14.   X = zeros ( size (x) );
  15.   
  16.   for( k in 1:x.n )
  17.   {
  18.     for( n in 1:x.n )
  19.     {
  20.       X[k] = X[k] + x[n]*exp( (-1i*2*pi*(k-1)*(n-1))/x.n );
  21.     }
  22.   }
  23.   
  24.   return X;
  25. };
  26.